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ULTRASONIC ELASTOGRAPHY WITH ANGULAR COMPOUNDING 

CROSS-REFERENCE TO RELATED APPLICATIONS 

[0001] This application claims the benefit of U.S. Provisional Application 
60/464,678 filed April 22, 2003 hereby incorporated by reference. 

STATEMENT REGARDING FEDERALLY 
SPONSORED RESEARCH OR DEVELOPMENT 

[0002] 

BACKGROUND OF THE INVENTION 
[0003] The present invention relates to ultrasonic medical imaging and in 
particular to an improved apparatus and method for making ultrasonic elastography 
images. 

[0004] Elastography is a new imaging modality that reveals the stiffness 
properties of tissue, for example, axial strain, lateral strain, Poisson's Ratio, Young's 
Modulus, and other common strain and strain-related measurements. The strain 
measurements may be collected over an area and compiled as a two-dimensional 
array of data, which may then be mapped to a gray or color scale to form a strain 
"image". Analogously, strain measurements may be collected over a volume 
displayed either three-dimensionally or as a series of stacked two-dimensional 
images. 

[0005] In quasi-static elastography, two images of the tissue may be obtained by 
the ultrasound device in two different states of compression, for example, no 
compression and a given positive compression. The tissue may be compressed by 
an external agency such as a probe or the like, or may be compressed by muscular 
action or the movement of adjacent organs. Strain may be deduced from these two 
images by computing gradients of the relative local shifts or displacements in the 
images along the compression axis. Quasi-static elastography is analogous to a 
physician's palpation of tissue in which the physician determines stiffness by 
pressing the tissue and detecting the amount the tissue yields under this pressure. 
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[0006] In dynamic elastography, a low-frequency vibration is applied to the 
tissue and the tissue vibrations are measured, for example, using Doppler detection. 
[0007] Ultrasonic elastographic images may have considerable image noise 
arising from a number of sources. In quasi-static elastography, one noise source is 
so-called "correlation artifacts" which relate to the amount of overlap in the 
windows of correlation as may be reduced improve resolution. 
[0008] A number of techniques have been used to minimize these image 
artifacts, including: global temporal stretching in which the post-compression echo 
signal is stretched by a known compression factor before estimating axial strain, 
wavelet de-noising where wavelet techniques are used to smooth the displacement 
estimates in the wavelet domain without losing edge information, and multi- 
compression averaging in which multiple readings are taken with smaller 
compression increments and averaged together. 

[0009] Typically elastography produces a strain measurement only along the 
axis of compression. However, lateral strain may be of value both in deducing 
qualities like Poisson's Ratio and in countering the effects of lateral motion in de- 
correlating the axial displacement of the tissue. In addition, shear strain images can 
also be obtained. 

[0010] Two principal methods have been used to obtain lateral strain. The first 
described in U.S. Patent 6,270,459 interpolates between successive axial rays or 
echo signals to provide a basis for horizontal displacement measurement using a 
correlation technique. A second technique makes the assumption that the tissue 
being measured is incompressible and simply deduces lateral strain from the axial 
strain. 

SUMMARY OF THE INVENTION 
[0011] The present invention provides a method and apparatus to significantly 
decrease image artifacts in elastographic imaging by using multiple angle 
acquisitions combined either before or after strain is extracted. The inventors have 
determined that even a relatively small angular difference between acquisitions 
provides the necessary statistical independence of the measurements to significantly 
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reduce image noise. This angular compounding concomitantly provides both axial 
and lateral strain measurements. 

[0012] Specifically, the present invention provides an ultrasonic elastography 
system including an ultrasonic transducer system adapted to provide a set of echo 
signals from at least two angles from a plurality of voxels in a region of interest. 
The set of echo signals includes first echo signals taken with the region of interest in 
a first compressive state (e.g. no compression) and a second echo signal taken with 
the tissue of the region of interest in a second compressive state. A processor 
combines at least portions of the set of echo signals from different angles to produce 
at least one compounded strain measurement for voxel in the region of interest. 
[0013] Thus is it one object of the invention to provide a simple method of 
decreasing the statistical fluctuation in elastographic measurements to improve 
image quality and quantitative accuracy. 

[0014] The processor may use the angle of the ultrasonic echo signal underlying 
the strain measurement to convert the strain measurement to an equivalent strain 
measurement along a predetermined strain axis and may combine the converted 
strain measurements to produce the compounded strain measurement. 
[0015] Thus is it another object of the invention to provide a method of 
combining strain measurements made at different angles to improve strain 
measurements made along a given angle by converting the multiple strain 
measurements to a common strain axis before combination. 

[0016] The predetermined strain axis may be the axis of compression or may be 
an axis lateral to that axis. 

[0017] Thus it is another object of the invention to provide strain measurements 
along arbitrary axes. 

[0018] The angle of compression of the tissue may remain constant for the 
different angles of echo signal. Alternatively, the angle of compression may change 
with changes in the angle of the echo signal. 

[0019] Thus it is another object of the invention to flexibly provide for single or 
multi-angle compressions to be used in extracting additional information from the 
tissue. 
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[0020] These particular objects and advantages may apply to only some 
embodiments falling within the claims and thus do not define the scope of the 
invention. 

BRIEF DESCRIPTION OF THE FIGURES 
[0021] Fig. 1 is a simplified block diagram of an ultrasound scanner suitable for 
use with the present invention; 

[0022] Fig. 2 is cross-sectional view of tissue showing collection of the 
ultrasonic echo signals obtained with the apparatus of Fig. 1 into echo signal sets 
having different ray angles; 

[0023] Fig. 3 is an alternative embodiment of the ultrasonic transducer system of 
Fig. 1 for multi-angled transducer and/or compression plate movement; 
[0024] Fig. 4 is a simplified representation of two echo signals obtained in 
different compression states showing their relative displacement; 
[0025] Fig. 5 is a signal flow diagram showing processing of the sets of Fig. 2 
into angularly compounded elastographic images showing axial and lateral strain; 
[0026] Fig. 6 is schematic representation of an individual voxel of tissue under 
compression showing angles used for the computation of the weighting factors of 
Fig. 5; 

[0027] Fig. 7 is a flow diagram similar to that of Fig. 5 showing expansion of 
the technique to multiple compression angles; 

[0028] Fig. 8 is a perspective representation of the acquisition of scan data over 
a three-dimensional region of interest for parametric imaging of a volume rather 
than a single plane; and 

[0029] Fig. 9 is an ultrasonic transducer showing scanning using a multi-element 
phased-array linear or curvilinear transducer. 

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENT 
[0030] Referring now to Fig. 1, an elastographic imaging system 10 suitable for 
use with the present invention may employ an ultrasonic imaging machine 1 1 alone 
or in combination with an external computer 30. Generally, the ultrasonic imaging 
machine 1 1 provides the necessary hardware and/or software to collect the series of 
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ultrasonic echo signals that can be further processed as will be described by a 
processor 33 held within the ultrasonic imaging machine 1 1 in the external computer 
30. 

[0031] An ultrasonic transducer 12, associated with ultrasonic imaging machine 
11, may transmit multiple ultrasonic beams 14a, 14b, and 14c (only three shown for 
clarity) toward a region of interest 18 within a patient. Ultrasonic beam 14b, which 
typifies the ultrasonic beams 14, produces echoes along different measurement rays 

16 acquired at different angles within each ultrasonic beam 14. Different rays 16 
will pass through different volume elements (voxels) 26 within the region of interest 
18 and generally a given voxel 26 will have multiple rays 16 passing through it at 
different angles from different ultrasonic beams 1 4. 

[0032] The echoes along each of the rays 16 will be received by the transducer 
12 and converted to electrical echo signals 17 acquired through the interface 
circuitry 22 of the ultrasonic imaging machine 1 1 . The interface circuitry 22 may 
provide amplification, digitization, and other signal processing of the echo signals 

1 7 as is understood in the art. These digitized echo signals may then be transmitted 
to a memory 35 for storage in subsequent processing by processor 33 as will be 
described. 

[0033] Ultimately, the processed echo signals 17 are used to construct an image 
displayed on a graphic display 32. Input commands from an operator may be 
received via a keyboard 34 and/or a cursor control device 36 such as a mouse 
attached to the processor 33 as is well understood in the art. 

[0034] In one embodiment, the ultrasonic imaging machine 1 1 may be an Aloka 
SSD 2000 scanner employing a phased array transducer with a center frequency of 
3.5 MHz with a 6 dB band of 640 percent. Digitized echo signals from this 
ultrasonic imaging machine 1 1 may be captured by a Gage Applied Science 12100 
A/D board and provided to the computer 30 for processing. 
[0035] In one embodiment, the ultrasound imaging machine 1 1 may be a 
scanner that provides an interface that enables the user to utilize beam-steering on a 
linear array transducer. The ultrasound beam can be steered at different angles as 
shown in Fig. 1 . The raw radio frequency data obtained at each beam-steered angle 
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is stored on the ultrasound system for either off-line or real-time processing on the 
ultrasound system itself. The strain images or elastograms will be displayed 
alongside the ultrasound B-mode sonograms or superimposed on the ultrasound 
sonograms. 

[0036] Generally, as shown in Figs. 1 and 9, the ultrasound transducer 12 may 
be a single element transducer manually steered to acquire echo signals along 
different rays 16 of each ultrasonic beam 14, or be a multi-element ultrasonic 
transducer 12 electronically steered by phased array operation to acquire echo 
signals along different rays 16 of each ultrasonic beam 14. The ultrasonic 
transducer 12 may be manually repositioned for each of the different ultrasonic 
beams 14 or automatically moved by means of a motorized stage 42 receiving 
position signals 41 from the processor 33 via the interface circuitry 22. As shown in 
Fig. 2, a multi-element linear or curvilinear array ultrasonic transducer 12 may be 
steered to collect a series of angulated beams with minimal or no physical 
translation. Mechanical translation and angulation may also be used with multi- 
element phased array transducers. 

[0037] A position sensor 19 may be attached to the ultrasonic transducer 12 to 
obtain position data 43 indicating the position and orientation of the ultrasonic 
transducer 12, or in the case where a motorized stage 42 is used, position data 43 
may be extracted from the position signals 41 as the motorized stage 42 is moved. 
[0038] Referring still to Fig. 1, during acquisition of the ultrasonic echo signals 
17, compression is applied to the patient via a compression plate 44 oriented to 
compress the patient along a transmission axis 46 generally aligned with the center 
rays 16 of the ultrasonic beams 14. The compression plate 44 has a central aperture 
allowing the ultrasonic transducer 12 to transmit ultrasound therethrough. The 
compression plate 44 may be manually operated or may be operated by means of a 
compression signal 47 obtained from the interface circuitry 22 and received by 
actuators on the compression plate 44. Alternatively compression may be applied by 
the ultrasonic transducer 12 itself or other forms of compression may be used as are 
understood in the art. Alternatively, compression can be applied using the 
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transducer itself, or with the attached compression plate in a freehand manner or 
using a motorized hand-held compressor for controlled compressions. 
[0039] Referring now to Fig. 3 in an alternative embodiment of the invention, 
the motorized stage 42 of Fig. 1 may be replaced with an arcuate gantry arm 50 
along which the ultrasonic transducer 12 may travel in an arc about the patient in 
region of interest 18, ideally over approximately 180 degrees. In this case ultrasonic 
beams 14 orbit about a center rather than translate as shown in Fig. 1. A 
compression plate 44 may be mounted with the ultrasonic transducer 12 to apply 
compressions along the transmission axis 46 or may be fixed at one angle 
independent from the transmission axis 46. Again, position signals 41 may control 
the position of the ultrasonic transducer 1 2 and compression signals may control the 
compression plate 44 and position data 43 provided to the processor 33. 
[0040] Referring now to Fig. 2, echo signals 17 along rays 16, from either of the 
embodiments of Figs. 1 or 3, may be optionally collected in projection sets 56 (only 
two projection sets 56a and 56b are shown for clarity) such that echo signals 17 in 
each projection set 56 are along parallel rays 16. Each projection set will obtain two 
different echo signals 17 and 17' (the latter indicated by dotted lines), one 
corresponding each to two different compression states of the compression plate 44. 
[0041] The angle between measurement rays 16 of each projection set 56 may 
differ by as little as 0.75 degrees and still provide sufficient independence of 
measurement to reduce the statistical deviation in the measured strain. Nevertheless, 
higher degrees of angular separation may also be used and in a preferred 
embodiment, for example, forty-five different measurement sets 56 may be obtained, 
each with one degree difference between them. Clearly higher angular differences 
between the rays 16 or the measurement sets 56 produce even greater statistical 
independence in the measurements such as will help reduce noise. The invention 
contemplates that angular separations of five degrees or greater may also be 
practical over an angular ranges of up to 180 degrees. An optimal angular 
separation exists that will be used when possible under practical imaging conditions 
using beam-steering on a linear array transducer. 
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[0042] Referring now to Fig. 4, generally within each projection set 56, an echo 
signal 17 acquired during a first degree of compression (e.g., no compression) will 
differ from an echo signal 17' acquired during the second degree of compression 
(e.g., a predetermined positive compression) by a dilation of the echo signal 17' 
along the horizontal axis or time axis with respect to echo signal 17. In the example 
of Fig. 4, the echo signals 17 and 17' are raw ultrasound data, however, similar 
analysis can be performed with B-mode or amplitude envelope data. 
[0043] As is understood in the art, for each pair of echo signals 17 and 17', 
relative displacement between the echo signals 17 and 17' at each point along the 
length of echo signals 17 and 17' are analyzed and through correlation or other 
techniques described generally in U.S. Patent Application 2003/30171672 Al, 
hereby incorporated by reference, to determine a displacement 60. This analysis is 
confined to a window 58 determined according to the desired resolution of strain 
measurement. This strain deduced from displacement 60 will be axial strain, that is 

r 

strain measured along the transmission axis 46 of the echo signal 17. As is generally 
understood in the art, each window 58 may be identified with different volume 
elements 26 through which the particular rays 16 of echo signals 17 and 17' pass. 
[0044] Referring now to Fig. 5 in a first step of the present invention, echo 
signal sets 62 are acquired, each including echo signals 17 and 17' at different 
compression amounts for one beam 14. Different echo signal sets 62 will hold data 
of different beams 14 taken at different locations and/or different angles by the 
apparatus described with respect to Figs. 1 or 4. The angle information of each echo 
signal 17 is preserved within the echo signal sets 62 as recorded by the processor 33 
during acquisition of this data. 

[0045] In a next step, the echo signal sets 62 are optionally provided to a 
rebinner 64 which separates them into parallel projection sets 56 described above, 
each again having echo signals 17 and 17' associated with two sets of compression, 
however, rebinned among echo signal sets 62 so that the echo signals of each 
projection set 56 are along parallel rays. The particular transmission axis 46 of the 
rays associated with each projection set 56 is determined and used later. Rebinning 
is necessary only when sector or phased array transducers are used. 
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[0046] A B-mode image (not shown) may be extracted from any or all of either 
the projection sets 56 or echo signal sets 62 to produce a reference image optionally 
to be projected together with the elastographic image on display 32. 
[0047] The projection sets 56 may then be provided to a strain extracting block 
66 which extracts strain information from the projection sets 56 by comparing the 
echo signals 17 and 17' as discussed above with respect to Fig. 4 to produce a series 
of elastographic images 68 corresponding generally to each of the projection sets 56. 
The elastographic images 68 will be composed of pixels 70, each representing a 
strain vector 72 having a magnitude proportional to the displacement 60 described 
above with respect to Fig. 4, and an angle determined from the transmission axis 46 
described above. 

[0048] The elastographic images 68 on a pixel by pixel basis are then weighted 
by weighting factors W and W to extract weighted elastographic image pairs 74 
comprised of an axial strain image 76 and a lateral strain image 76'. In this 
extraction, the pixels 70 of the axial image 76 provide strain extracted to a 
predetermined axial direction 51 and the pixel 70 of the image 76' provides 
extraction along a lateral direction 5 1 ' perpendicular to the axial direction 51. It will 
be understood that axial direction 51 and lateral direction 51' may be arbitrarily 
selected, however, typically the axial direction 5 1 will be along the direction of 
compression if there is one compression direction, and for the embodiment of Fig. 3 
where the direction of compression is changing to a midrange of the gantry arm 50. 
[0049] Referring now to Fig. 6, generally the weighting factor used for this 
extraction will be a function of the angle 0 between the transmission axis 46 and the 
axial direction 5 1 . The weighting factor may be derived using a variety of methods, 
however, in the preferred embodiment, it is deduced by assuming that the voxel 26 
under compression along transmission axis 46 is an incompressible material, in other 
words, a material having a Poisson Ratio of approximately 0.5. In this case, a strain 
along transmission axis 46 may be multiplied by weight w(9) to produce the 

equivalent expected strain along transmission axis 46 according to the following 
equation: 
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w(0) = 



1 



(1) 



(cos 2 <9-vsin 2 6) 



[0050] where v is the Poisson's ratio is equal to 0.495 for incompressible media. 
Similarly, strain along transmission axis 46 may be multiplied by a slight variant of 
this weight W* to produce strain along lateral direction 51', per the following 
equation: 



[0051] Referring again to Fig. 5, the weighted strain vectors 72 of the pixels 70 
in each of the images 68 are thus resolved to axial strain vectors 73 and lateral strain 
vectors 73' in elastographic image pairs 74 allowing like-oriented strain vectors of 
each of the elastographic image pairs 74 to be to be averaged per process block 80 to 
produce two images, an axial strain image 82 and a lateral strain image 82'. This 
summing process reduces the noise in each of the images 82 by an amount roughly 

equal to where N is the number of images 68. 

[0052] The axial strain image 82 and a lateral strain image 82' may be mapped 
by means of mapper 84 to the graphic display 32. For example, shades of gray may 
be used to map magnitudes of the axial strain and different colors may be used to 
map magnitudes of the lateral strain. Alternatively, the axial displacement and strain 
image 82 and lateral displacement and strain image 82' may be combined to produce 
shear strain images. Improved reconstruction of Young's Modulus images can also 
be obtained with more accurate normal and shear strain tensor information. A B- 
mode image may also be displayed next to or superimposed over the lateral image or 
switched with the lateral and axial images to provide a point of reference for the 
particular strain values. 

[0053] Referring now to Fig. 7, the weighting block 83 may provide a more 
complex weighting system when the axis of compression may move. In this 
embodiment, weights from weighting block 83 are applied to the images 68 to 
produce a large number of compression identified elastographic image pairs 74 each 
associated with a different angle of compression axis. These may be combined per 
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process block 80 as before to produce an axial strain image 82 and a lateral strain 
image 82' for each compression axis. Each of the axial strain image 82 and lateral 
strain image 82' may, in its own right, produce an image for a particular 
transmission axis 46 or they may be combined by a combiner 90 to produce some 
directionless quantity such as Poisson's Ratio or Young's Modulus from each of 
these different measurements, such as made to be displayed on graphic display 32. 
[0054] Referring now to Fig. 8 for reasons of clarity, the invention has been 
described with respect to voxels 26 aligned in a single plane corresponding to a 
plane of the image produced on graphic display 32. However, it will be understood 
the essential principle of summing together the echo signals taken at different angles 
to improve elastography measurements may occur by moving the ultrasonic 
transducer 12 so as to collect multiple ultrasonic beams that differ not only by their 
angle within the plane but also in angles over a three-dimensional curved surface so 
as to produce a volumetric image data that may be displayed one slice at a time or 
rendered as a three-dimensional object. 

[0055] It will be understood generally that the order of combining the images 
and extracting the elastic strain information is arbitrary and that in certain instances 
the images may be combined first and then strain data extracted or strain data 
extracted and then combined in an angularly sensitive way. 

[0056] It is specifically intended that the present invention not be limited to the 
embodiments and illustrations contained herein, but include modified forms of those 
embodiments including portions of the embodiments and combinations of elements 
of different embodiments as come within the scope of the following claims. 
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